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£C) , Abstract 

■ A calculation of n° — > e + e~7 in next-to-leading order in chiral per- 

fN| ' turbation theory of two flavour case extended by virtual photons is 

presented. The whole kinematic sector is covered and realistic experi- 
' mental situation is discussed. 

P-i. 1 Introduction 

P-i 

With a branching ratio of (1.198 ± 0.032)% 0] the three body decay vr° -> 
e + e~7, first considered by R. H. Dalitz more than fifty years ago 0, is 
the second most important decay channel of the pion. The dominant decay 
^ . mode, 7T° — > 77, with its overwhelming branching ratio of BR = (98.798 ± 

0.032)%, is deeply connected with this three body decay. The other decay 
channels connected with the anomalous 7r°77 vertex, like 7r° — > e + e~ and 
7T° — > 2(e + e~), are suppressed approximately by respective factors of 10 -7 
and of 10 -5 . As already alluded to, the hadronic aspects of the Dalitz decay 
of the neutral pion are related to the axial anomaly. Actually, the main 
features of this process can be discussed within the framework of Chiral 
Perturbation Theory (ChPT - for reviews see @). In the case at hand, we 
need only to consider the case of two light quark flavours, u and d. We 
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Figure 1: Kinematics lay for 7r° — ► 77 and 7r° — > e + e 7 

shall however use the extension of ChPT to electromagnetic interaction, i.e. 
including effects of virtual photons ||], treating e.g. the pion mass difference 
as a p 2 effect at lowest order. 

2 Kinematics 

One may separate the contributions to the ir° — > e + e~7 process into two 
main classes. The first corresponds to the Feynman graphs where the 
electron-positron pair is produced by a single photon (Dalitz pair). The 
leading contribution, of order a 3 , to the decay rate belongs to this one- 
photon reducible class, which involves a semi-off-shell 7r°77 vertex. The 
second class of contributions corresponds to one-photon irreducible topolo- 
gies. They start with the radiative corrections to the 7r° — > e + e~ process, 
which involves the doubly off-shell 7r°77 vertex. Their contributions to the 
decay rate are thus suppressed, being at least of order a 5 . For the time 
being, we shall therefore not take them into account. A complete discus- 
sion will appear elsewhere H. The general expression of the one-photon 
reducible part of the ir°(p) — > e + (qi)e~ ((72) j(k) amplitude is of the form 

M ljR (n° - e+e" 7) = M v <f"> ( ~* } - - 1 - 

Wi+92) 2 l + n((gi +q 2 ) 2 ) 

x u(q2)(ie)T p (q 2 , q^v^) . (1) 

In this expression, M^p^k) represents the semi-off-shell 7r°77 amplitude, 
defined as 

M^(p,k) = ( 7 (A;);out|^(0)|7r (p);in) = A((p - k) 2 )e v ^e u (k)*k a pP , (2) 

where the second expression follows from Lorentz invariance and parity 
conservation, while denotes the electromagnetic current. Furthermore, 
n ((<7i + <Z2) 2 ) denotes the renormalized vacuum polarization function, with 
n(0) = 0, and Tp(q 2 , qi) stands for the vertex function, renormalized in the 
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on-shell scheme, of the current j p . Assuming again parity conservation, one 
may express T p (q 2 ,qi) in terms of two form factors 

Tpfe, 9i ) = F 1 (( gi + q 2 f) lp + F 2 (( gi + q 2 f) ^ P Mi + Q2) '■ (3) 

The Dirac form factor i<i is thus normalized by -Fi(O) = 1, while the Pauli 
form factor gives the anomalous magnetic moment of the electron, ^(0) = 
a e . Finally, since (qi + q2) fl -M-fi(p, k) = 0, the photon propagator reduces to 
the contribution of the g up piece. 

It is convenient to introduce the following kinematic parameters 



(gi + Q2) 2 
2 ' 

2p • (gi - q 2 ) 
m^ (l - x) 



u z < x < 1, 



-P<y<P, 



4m 

2 ' 
m 



/9 



1 . 



(4) 
(5) 
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One may then write the partial decay rate, normalized to the tt° 
decay rate To = (l/647r).A(0)m^o , for the one photon reducible contribution 
in terms of these parameters as 



1 iff _a A(m 2 n0 x) 1 
T^dxdy ~ 7T A(0) l + n(m 2 x) 



2(1 - x) 3 r, ,2 , 2 

-||Fi(m; z)| (z/ + xy 2 + x) 



4x 2 



+ I F 2 (m» 1 2 ^ J 2 + X) + 4Re(F 1 (m^)F 2 * (m^ x))x} , 



(6) 



or, by integrating out the parameter y 



1 dT a 
Tq dx 7T 



1 



+ \F 2 (mnx] 



A(0) 1 + n(m 2 ; 
,2 x(2^ 2 + x) 



'jl-xf 
3x 2 



/?{|^iKo 



x)| 2 (z^ 2 + 2x) 



+ 6Re(Fi(m 2 x)F 2 *(m 2 x))x}. 



(7) 



At lowest order in the fine structure constant we have II = 0, F\ = 1 and 
F 2 = 0. We shall discuss A and the hadronic part of II in the next section. 

3 Parameterization of Hadronic Part 



We are going to study the problem in two flavour Chiral Perturbation Theory 
extended by virtual photons H,f?j up to one loop, i.e. for this type of 
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Figure 2: Leading order and pion loop contributions to the amplitude A4i^r oc 
.4(1 — Tlhad)- To simplify the figure, we have omitted the e + e _ 7-vertex. The 
graphs in the second row correspond to the corrections to A, while the third row 
contributes by the factor 1 — Tlhad- 

process up to the 0(p e ) where p can represent momenta relevant for the 
decay or electromagnetic constant e. The lowest order effective Lagrangian 
can be written as 

£ (2) = ^l^u+^u + x +u + u+ x) _ If^f^ + ZF\Q R UQ L U+) (8) 



with 



i± ( vr° y/2: 



7T 



+ 



Further chiral invariant Lagrangians of higher order (0(j> 4 ) and 0(p e )), so- 
called Wess-Zumino anomaly Lagrangian (of 0(p 4 ) order) and also the other 
details of calculation within ChPT can be found in above references. 

We will calculate the hadronic part of one photon reducible amplitude of 
the decay tt° — ► e + e~7, as it is depicted in Fig.0, i.e. the value of .4.(1— tlhad) 



in the terms of the previous section. With definition a = J 1 — 4m? ± /I 2 one 
can obtain from the one photon irreducible graphs in Fig.0 



A(l 2 ) 



N c e 2 r I 2 1 - - 2 



12vr 2 F^o I 6F 2 (4vr) 2 



log^ + i-2a 2 

m „ 



ex 3 log - — - - Aol + 5a } , (10) 
cr + 1 J J 
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while the one particle reducible graphs contribute by the factor (1 — Tlhad), 
where the hadronic contribution to renormalized vacuum polarization func- 
tion is 



487T 2 



| + 2c7 2 + ( j 3 log^— \ 
° a + 1 



(11) 



Here we have used standard handling with infinities (A) for 0(p e ) constants: 

(12) 



Aa = Aa r + AaX 
5a = 5a r 



with renormalized quantities to be finite and 

,2 



Aa 



Aa 
32^2 



, , mz.o \ 

(Aa + log^j-J, Aa 



12vr 2 ' 



(13) 



The 0(p e ) constants used in this text could be connected with form of 
and U respectively as 



Aa = -F$(A 2 - 4A 3 ) = -F 2 c 13 



A 32 2 

= T 71 " 

64 2 

= 7T 

3 
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(A 2 - 2A 3 - AA A )m^ Q + — {A 4 + 2A 6 )2rhB 



(en - 4c 3 - 4c 7 )m 2 + -(5c3 + c 7 + 2c 8 )2mB 



(14) 



(15) 



The influence of virtual photon is covered in ( jig ) in the neutral pion 
decay constant F n o M. Unfortunately this constant is not known very ac- 
curately F n o = 92 ± 4 Q and so in practice the charged pion decay constant 
is used F v = 92.4 



2 

m + 



2 

m + 



F n o = F n (1 - ^-f^ log -f± - ^ [(3 + IZ)h -fZk 2 ~ 3k 3 ~ ±Zh 

(16) 

To cover this section we will introduce a slope parameter a which is a 
parameter of linear expansion of A(l — Tlhad) in x: 



A(l - U had ) ~ (1 + ax) 



(17) 



4 Calculation at e 3 (l + p 2 + e 2 ) order 

Now, to the above hadronic contribution we have to add also the QED cor- 
rections, which are in our case up to one loop level: lepton part of vacuum 
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polarization function, vertex, bremsstrahlung and 2-photon triangle and box 
corrections. In the first approximation we will neglect 2-photon loop dia- 
grams (for a polemic whether this is justified see also [|K]). For one-photon 
exchange the following interesting relation holds (c.f. |pA|) 

r(vr° -» e+e- 7 ) = 2 / — -Imn(s)r (s). (18) 
Jo x TT 

with the lepton part of spectral function of the photon propagator II(s). 
With 0(p e ) constants set equal zero and therefore (from (|To|) , dTI| ) and 
(|T7D) a = —0.0024 we can obtain numerically 

pfvr — ► e + e~^) 
K — „ -^ = 0.01185 + 0.00010 = 0.01196, (19) 

r(7r u — > 77) 

where the second number represents loop-correction contributions. Com- 
paring this with PDG's number T ex P(7r -> e + e-7)/T(vr -> 77) = (1.213 ± 
0.033)% one can immediately see that the total decay rate is not very suit- 
able for probing the smooth effects of the higher orders (e.g. tuning ChPT 
constants) . 

Therefore we will turn our attention to the differential decay rate, which 
can be calculated in one-photon case using the relation (obtained from (|l8|)) 

: -n a (x)r (x) + / —U b (E,x)T (E), (20) 



dx 



where we have divided the spectral function to the virtual and real part 

-ImU^(p)=U a (p 2 ) + U b (p 2 ). (21) 

7T 

From the first plot in Fig.||| we see that while the radiative effects for 
the total decay rate is small, the situation in the differential decay rate is 
opposite. 

However, this approach is useful for a situation when the detector is 
entirely blind to photons. In a more realistic situation, the experimental 
resolution is such that it is possible to detect photons if only above the 
(small) energy AE. For this case we define 5(x,y) and S(x) 

1 d 2 r rad e , 1 d 2 r 1 dr rad r/ i dr 

=S(x)^r^-. (22) 



To dxdy Tq dxdy Tq dx Tq dx 

In Figj^ (the second plot) we depicted how the 5(x) changes with parameter 
aE. 
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aE = WMeV 

Figure 3: Differential decay rate j^^jj— - (from (|(])) and the dependence 
of 5(x) on aE 

Up to now we have systematically neglected the graphs which arise from 
the two-virtual-photons exchange. We have performed this calculation and 
we found interesting behaviour of contribution to S(x), particularly when x 
is close to 1. In this domain 

<5(x~l,y)«0.1, (23) 

i.e. it represents approximately 10% of leading order and from the previous 
we could conclude that it is a important contribution to differential decay 
rate, which cannot be neglected. However, it will be difficult to test this 
area experimentally. Actually for 

*>!-— (24) 

the photon is effectively immeasurable by definition and precisely this is the 
domain where the interesting situation of ( |2~3| ) occurs. Let us also stress that 
in this part of calculation we have neglected the electron mass and thus, at 
this level, it is not possible to compare and use in our consideration process 
7T° — > e + e~ which is proportional to m e (for more details see 

5 Conclusion 

We have discussed the Dalitz decay ir° — ► e + e~7 in two flavour ChPT and 
set the hadronic part within this effective theory. Using PDG's value for the 
slope parameter (see (|l7|)) a = 0.032 ± 0.004 [|| one can estimate following 
linear combination of 0(p & ) constants 

m 2 1 

— ^ 7- NO Aa -a5a = a + 0.0024 
6F 2 (47r) 2 
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(cf. ([D]) and (|l4|)). We approve that the effects of higher orders in total 
decay rate are smaller than the present experimental error (|l9|). However, 
differential decay rate is much sensitive to these effects. 
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